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Abstract 

The Lorentz covariant classical and quantum statistical mechanics and thermodynamics of an 
ideal relativistic gas of bradyons (particles slower than light), luxons (particles moving with the 
speed of light) and tachyons (hypothetical particles faster than light) is discussed. The Lorentz 
covariant formulation is based on the preferred frame approach which among others enables consis- 
tent, free of paradoxes description of tachyons. The thermodynamic functions within the covariant 
approach are obtained both in classical and quantum case. 
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I. INTRODUCTION 



Ideal gases are one of the most important model systems of the nonrelativistic statistical 
mechanics and thermodynamics. Examples range from equations of state for classical gases 
to description of electrons in metals and superconductors. In spite of the fact that the studies 
of a relativistic gas of massive particles (bradyons, also called tardyons) go back to 1911 l|, 
the relativistic statistical mechanics and thermodynamics are far from complete. On the 
one hand, the reason are apparent limited applications of the relativistic thermodynamics. 
For example, in opinion of Ter Haar and Wergeland [2]: "At extremely high temperatures 
relativistic effects may of course be important. Then, however, matter behaves as mixture 
of ideal gases and this limiting case poses no problem. By and large, a relativistic theory 
of heat seems, therefore, to be of little practical importance." Nevertheless, the arguments 
of Ter Haar and Wergeland evidently fail for luxons and tachyons which are relativistic 
particles regardless of the concrete value of the temperature. Furthermore, as pointed out 
by Aragao de Carvalho and Goulart Rosa [3], a fully relativistic treatment is required by 
astrophysical systems such as white dwarfs and neutron stars. On the other hand, the 
development of the relativistic statistical mechanics and thermodynamics was slowed down 
by the lack of the covariant formulation. In particular, we point out different transformation 
rules of the relativistic temperature suggested by Einstein, Planck and von Laue, by Ott 
and by Landsberg . The only formulation working in the case of tachyons, based on the 
nonstandard (absolute) synchronization scheme, was introduced in a very recent paper J]. 
In this work we study the Lorentz covariant statistical mechanics and thermodynamics of the 
relativistic ideal gas of bradyons, luxons and tachyons. In section 2 we recall the formulation 
of special relativity based on the absolute synchronization. Section 3 is devoted to the 
classical relativistic ideal gas. In particular, we derive the covariant form of thermodynamic 
functions. In section 4 we discuss the quantum statistical mechanics and thermodynamics 
of the relativistic ideal gas. Besides derivation of the covariant forms of thermodynamic 
functions we also discuss the classical limit. 
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II. ABSOLUTE SYNCHRONIZATION SCHEME 



In this section we recall the basic facts about the formulation of special relativity with the 
help of the absolute synchronization of clocks [5]. Among others this approach provides a 
correct description of tachyons. Tachyons are hypothetical faster than light particles. Besides 
their intriguing theoretically predicted properties |6j, tachyons take attention of physicists 
because they are candidates for the dark matter [7] and dark energy [8j]. Moreover they 
appear in brane theories such as the brane excitations as well as in cosmological models (so 
called rolling tachyon models) [9j. For this reason it is interesting to investigate statistical 
and thermodynamical properties of tachyonic gas. Unfortunately, standard description of 
tachyons is plagued by number of inconsistencies. Typical examples of such difficulties are 
the causal paradoxes (tachyon anti-telephone jlo|). the problem of so called transcendental 
tachyon (the space of tachyon velocities is not a Lorentz group carrier space 0, [ll|) and 



vacuum instability on the quantum level 



12|. 



As was stated many years ago by Sudarshan Uj a consistent description of tachyons 
demands a preferred reference frame on the fundamental level. However, this means that 
the relativity principle is necessarily broken in such a case. This causes an apparent conflict 
with the standard Lorentz group transformations in the Minkowski space-time. To overcome 
this difficulty let us notice that introduction of the inertial preferred frame (PF) means that 
we should realize the Lorentz group not only on the space-time coordinates but also on the 
four-velocity of the PF as seen by inertial observers. This gives us the necessary freedom 
to reconcile breaking of the relativity principle and simultaneously to preserve Lorentz co- 
variance. Such a realization of the Lorentz group was given in [5| and it has an elegant 
explanation in terms of the bundle of frames as we 1 as ;he physical interpretation in terms 
of the absolute synchronization scheme for clocks (5,(1^, 3- m particular in \^ a consistent 
classical and quantum description of tachyons was built in this framework, without of the 
above mentioned inconsistencies. It is important to stress that for massless and massive sub- 



luminal particles (luxons and bradyons respectively) this scheme is completely equivalent to 
the Einstein synchronization scheme (so called convention of synchronization [5]) whereas it 
provides a consistent description of tachyons. As an important application of the absolute 
synchronization in quantum mechanics the covariant relativistic position operator was re- 
cently introduced 14]. In the PF this operator reduces to the well-known Newton- Wigner 



3 



operator. Finally, let us notice that the absolute synchronization is the most natural in 
cosmology because our universe distinguishes a frame (cosmic background radiation frame) 
as well as is flat on the large scales. Let us also recall that some recent theoretical inves- 
tigations of quantum gravity and extremely high energy phenomena predicts existence of a 
preferred frame. For example the loop gravity describing quantum gravitational phenom- 
ena at the Planck scale predicts the existence of a preferred frame and even incorporates 
quantum scale in the corresponding effective Lorentz group realization (the so called DSR 



theories 



15| and the Einstein-aether theories 16]). Furthermore, in approach by Kostelecky 



a breaking of Lorentz symmetry is assumed via specific field interactions [17| • Consequently, 
in this approach there exists effectively a preferred frame of reference too. 

Let us recall briefly the main results related to the description of tachyons in the frame- 
work of the absolute synchronization. As was stated above, in the absolute synchronization 
scheme the Lorentz transformations are realized simultaneously on the both coordinates in 
the inertial reference frames and velocity of PF, namely the contravariant transformation 
rules between the frame O u > and O u are of the form 

x'(v!) = D(A,u)x(u), (2.1) 
u' = D(A,u)u. (2.2) 

where A is an element of the Lorentz group, u is the four-velocity of the PF with respect to 
the inertial observer. In (2.1) rotations are realized standardly i.e. D(R,u) = R, R e SO(3), 
while boosts are w-dependent: 



D(w, u) 



w 







U°W Cg> U T 



(2.3) 



where is four- velocity of the primed frame with respect to the unprimed one, and w ® u T 
designates the direct (Kronecker) product of the column vector w and the row vector u T . 
Hereafter the three-vector part of the covariant (contravariant) four-vector (a M ) will be 
designated by a (a). It is very important that the matrices D(w,u) are block-triangular, so 
the time coordinate is rescaled only by a positive factor. Namely, 

x '0 = J_ x o, (2.4) 
This enables us to avoid all inconsistencies which plague the standard approach to tachyons 

n n n 

[6|, HO, LL3J because the notion of the instant-time hyperplane is invariant under Lorentz 



transformations. On the other hand, the form of (2.4) allows to solve problems arising in 
the standard approach, of covariance of some relativistic observables such as the position 
operator mentioned earlier, related with mixing the time coordinate with the spatial ones 
in the transformation rule for x°. From the technical point of view the fact that in the 
standard Einstein synchronization the time coordinate is mixed with the spatial ones is 
connected with the non-triangular form of the Lorentz matrix given by (compare with the 
eq. (2.3)) 



A(io E ) 



T 

-w E 



1 + 



(2.5) 



where the subscript E designates the Einstein synchronization coordinates. It is easy to see 
that the relationship between the coordinates in the Einstein and the absolute synchroniza- 
tion is of the form 

x E = X° + U°U-X, X E = X, 



u E = -o, u E = u. 



Therefore the difference between both synchronizations lies in the definition of the time 
coordinate. Notice that the time lapse in a fixed point x is the same in both synchronizations 
(i.e. if dx = 0, then dx E = dx°). Moreover, the transformation (2.1), similarly as the 
standard Lorentz transformation preserves the notion of an inertial frame. The light velocity 
over a closed path is frame-independent as well. Furthermore, from (2.6) it follows that 
both schemes (Einstein and absolute) are equivalent for velocities less or equal to the light 
velocity (i.e. for bradyons and luxons), but for superluminal velocities (i.e. for tachyons) 
this equivalence is broken. From the technical point of view such nonequivalence is related 
to a singularity of the relationship between tachyon velocities (derived from (2.6)) in both 
synchronizations. 

The invariant Minkowski line element in the new coordinates in the frame O u has the 
following form: 

ds 2 = gflu (u)dx»dx\ (2.6) 

where the (covariant) metric tensor is frame-dependent (but not point-dependent!) and is 
given by 

1 



u°u 



U°U T 



-I + (u°) 2 u(g)u 



(2.7) 
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The contravariant metric tensor is of the form 



br(«)] 



(u ) 2 


U°U T 


u°u 


-I 



(2.8) 



Notice that from (2.9) it follows that the space metric is Euclidean in each frame O u i.e. 
dl 2 = dx 2 . Furthermore, u^u^ = 1, and space part of the covariant four-velocity is equal to 
zero i.e. u = in each frame, consequently u° — 1/uq ■ One can also easily check that we 
have 

W " ul = L (2 ' 9) 

It is also useful to express the four- velocity w in terms of the velocity v = w/wq of the 
primed frame with respect to the unprimed one by means of the formula 

w° = ; 1 (2.10) 



1 + u^u-v) — v z 



obtained from the relation w^w^ = 1. 



The following remarks are in order. The dispersion relations for four-momentum p M are 
given by 

p^Pn = m 2 , (bradyons) (2-H) 

p> M = 0, (luxons) (2.12) 

p^p^ = —m 2 . (tachyons) (2-13) 

We can solve these equations with respect to p°. It follows that 



p° = m°y (u-p) 2 + p 2 + m 2 , (bradyons) (2-14) 

p° = u°J(u-p) 2 +p 2 , (luxons) (2.15) 



p° = u° J (u-p) 2 + p 2 — m 2 . (tachyons) (2-16) 



We point out that p is in general different from the energy p given by 

p = —q[— u-p + yj (u-p) 2 + p 2 + m 2 ], (bradyons) (2-17) 

Po = ^ [—u-p + \J (u-p) 2 + p 2 ] , (luxons) (2.18) 

Po = u-p + yj (u-p) 2 + p 2 — m 2 ]. (tachyons) (2-19) 
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Now, the volume element dV = d 3 x = dx l A dx 2 A dx 3 transforms (under condition that 
dx° = 0) according to 

dV = w°dV, (2.20) 

so 

V = w°V, (2.21) 

where V = J. = . dV. Similarly, taking into account that the Lorentz invariant momentum 
measure has the form 

dfi(p) = 6{p°)5(p 2 - m 2 )d A p, (2.22) 

that is 

d 3 p 

dM = (2.23) 

where p is covariant momentum three vector, we deduce that d 3 p = dp±Adp2Adp^ transforms 
as d 3 p' = (l/w°)d 3 p. Therefore we obtain the following formula on the Lorentz invariant 
phase space measure: 

dT = -d 3 xd 3 p. (2.24) 

To complete the transformation rules discussed in this section we finally write down the 
following relations: 

P°' = ^P , u°' = — Q u°, u ' = w°u , (2.25) 
which are also used in the next sections. 



III. CLASSICAL RELATIVISTIC IDEAL GAS 

We now discuss the basic properties of the classical ideal gas in the absolute synchroniza- 
tion. Taking into account the transformation properties discussed in the previous section it 
is easy to deduce that the Lorentz invariant partition function for each particle is given by 

Q 

Z x = V J rf 3 pexp(-M^°), (3.1) 

where V is the volume of the system, (3 = l/(kT) and temperature transforms as T' = w°T 
under Lorentz transformations. It should be noted that in the preferred frame specified 
by Mo — 1 an d u = 0, Z x has the standard form. Furthermore, by means of the eq. (3.1) 
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14 



and the standard definitions of the thermodynamical functions it is easy to show [4| that 
temperature, internal energy, enthalpy, Helmholtz free energy and Gibbs free energy of the 
ideal gas transforms analogously to the volume (eq.(2.22)), whereas pressure, entropy and 
partition function are Lorentz invariant. Consider the case of the relativistic ideal gas of N 
noninteracting bradyons. In the case of bradyons the partition function (3.1) takes the form 

Z lb = V J d 3 pexp (- M o/3^/(-u-p) 2 + p 2 + m 2 ) , (3.2) 

where m is the rest mass of a particle and we set c = 1. Using the identities (A.l) and (A. 2) 
we find 

2 

Tfl 

Z lb = AnV^-K 2 (u Pm), (3.3) 

where K2(x) is the modified Bessel function (Macdonald function). In the preferred frame, 
when uq = 1, we obtain from (3.3) the well-known Jiittner result [lj. Furthermore, the 
partition function for luxons can be written as 

Z u = V j d 3 pexp (- Mo /5y / (n-p) 2 + p 2 ) , (3.4) 

which leads with the use of (2.10) to 

Finally, in the case of tachyons the partition function is given by 

Z lt = V J d 3 pexp (-u f3^(u-p) 2 + p 2 - m 2 ) . (3.6) 
Taking into account (A. 5) and (2.10) we get 

o 

771 

Z u = 47rV r ^-S , o, 2 («o/3m), (3.7) 
where Sq^Oe) is the Lommel function (see Appendix). In the preferred frame when uq = 1, 



iQ. 



the formula (3.7) reduces to that originally obtained by Mrowczyriski 

Now the partition function Z for ideal gas of N noninteracting particles is 

Z = ±Z». (3.8) 

The knowledge of the partition function enables calculation of thermodynamical quantities. 
The average energy U is related to the partition function Z by 

U = -p»Z- (3.9) 



Hence, using (3.8), (3.3) and (A. 2) we find the following formula on the average energy of 
the bradyon ideal gas: 

Ub ^JL +N ^kf^ . (3 . 10 ) 

p K 2 {u pm) 

We remark that (3.9) implies the transformation rule for the average energy of the form 
U' = w°U. From (3.10), (A. 2) and (A. 4) we find the following approximation of £/& for 

uof3m 1 

3N 

U b = -— + Nu m, uo/3m»l. (3.11) 
z p 

From (3.11) it follows that [/& approaches Nu^m for large (3 (see Fig. 1). We point out that 
in the second extreme case ito/3m <C 1, we have for an arbitrary thermodynamical quantity 
X 

lim X b)t = X h (3.12) 

following from the well-known fact that in the high-energy limit or equivalently high- 
temperature limit bradyons and tachyons behave as luxons. 

We now return to (3.9). An immediate consequence of (3.8), (3.9) and (3.5) is the 
following expression for the energy of the luxon ideal gas: 

U, - f . (3.13, 

Notice that U\ approaches zero as f3 approaches infinity (see Fig. 1). Referring to (3.12) it is 
also clear that (3.13) is the limit of (3.10) for uoj3m <C 1. Finally, taking into account (3.8), 
(3.9), (3.7) and (A. 7) we find that the average energy of the tachyon ideal gas is given by 

rr _ 3N AT UomS- ltl {uopm) 

U t — — o iV 5 7 n \ • (A i4 J 

P S 0i2 {u /3m) 
The relation (3.14) can be written in an equivalent form 

U t = -— + 3N — — - — — — 3.15 

P S Qt2 {uoPm) 

following directly from the second equation of (A. 8). In the preferred frame, when uq — 1, 
(3.15) reduces to the formula originally derived by Mrowczyhski [ijj]. The dependence of 
the average energy of the ideal relativistic gas on the temperature is shown in Fig. 1 and 
Fig. 2. On using (3.14), (A. 8) and (A. 10) we arrive at the approximate relation: 

2N 

Ut = —, iW3m»l. (3.16) 




FIG. 1: The plot of the average energy U per particle of the ideal gas of bradyons (dotted line), 
luxons (dash line) and tachyons (solid line) given by equations (3.10), (3.13) and (3.14), respectively, 
where no = 1, and m = 1 in (3.10) and (3.14), versus kT. The behavior of average energy near 
the absolute zero. 

Therefore, as with luxons, for large (3 the average energy U t approaches zero (see Fig. 1). 

As an immediate application of the obtained relations (3.10), (3.13) and (3.14) we now 
write down the formulas on the specific heat per particle C v related to the average energy 
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FIG. 2: Average energy U per particle of the ideal gas of bradyons (dotted line), luxons (dash line) 
and tachyons (solid line) given by equations (3.10), (3.13) and (3.14), respectively, where «o = 1 ; 
and m = 1 in (3.10) and (3.14), as a function of kT on a logarithmic scale. The parameters uo and 
m are the same as in Fig. 1. 
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Namely, we have for bradyons, luxons and tachyons, respectively 

uo/3mKi(uol3m) ( uo/3mKi(uof3m) 



a 



ub 



3 + (uopm) 2 - 3- 



C v t — k 



3 + 



K 2 {u Q f3m) 
C v i = 3/c, 

uoflmS^i t i(uol3m) — (« /5m) 2 S- 2i o(uof3m 
S 0i2 (u Pm) 



K 2 (u (3m) 



S 0j2 (u f3m) 



(3.18) 
(3.19) 

(3.20) 



The formula (3.20) can be written in an equivalent form 



a 



vt 



-k 



1 + 



3{uoPmS-i^(u f3m) — 5(u flm) 2 S_ 2ti (u f3m)} f 2,uol3mS-i^(uof3m 



S 0i2 (u l3m) 



+ 



S , 2 (u f3m) 



(3.21) 



following directly from (A. 8). Up to some typo in [18( (3.21) coincides in the preferred 
frame, when uq = 1, with the formula originally obtained by Mrowczyhski. We now return to 
(3.20). Taking into account (A. 8) and (A. 10) we obtain from (3.20) the following asymptotic 
relation: 



a 



vt 



In the case of bradyons we find 



C,, 



2k, 



2 ' 



u f3m ^> 1. 



Uo/3m ^> 1. 



(3.22) 



(3.23) 



We remark that the asymptotic relations (3.22) and (3.23) can be formally obtained from 
(3.17), (3.16) and (3.11). Furthermore, it follows from numerical calculations that, in oppo- 
sition to the case of bradyons when C v b is a decreasing function, the specific heat for tachyons 
has maximum (see Fig. 3). The occurence of the maximum was treated by Mrowczyhski [3] 
as a formal property of C v t- On the other hand, the maximum of the specific heat C v can be 



connected with the so called Schottky anomaly in a two-state system [19|. In the preferred 
frame, when uq = 1, the mass m in the argument Uofim of the function C v t is a counterpart 
of the energy E of a two-state system exhibiting Schottky anomaly, with a ground state of 
energy and excited state of energy E. Although a physical interpretation of the maximum 
of the specific heat C v t is not clear, nevertheless, the experience with the Schottky anomaly 
shows that the knowledge of the abscissa of the maximum of the specific heat would enable 
determination of the tachyon mass m. 
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FIG. 3: The plot of the specific heat of the ideal gas of bradyons (dotted line), luxons (dash line) 
and tachyons (solid line) given by (3.18), (3.19) and (3.20), respectively, where uq = 1 and k = 1, 

versus T/m 

We now discuss the entropy of the relativistic ideal gas. The Lorentz invariant entropy 
S can be expressed in terms of the partition function and the average energy. Namely, 

S = k\nZ + ^. (3.24) 

Using (3.24), (3.8) and the Stirling formula 

N\ « f^pj , (3.25) 
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S i = kN \ ln TJ^ + 4) (3-28) 



which holds for large N, we get 

S = kN\nZ 1 - N\nN + N + (3.26) 
From (3.26), (3.3) and (3.10) we obtain in the case of bradyons 

S b = kN ( In ^ + In K ^ U f m) + M °y i( I°! m) + A ■ (3.27) 

The formula on the entropy for the ideal gas of luxons such that 

8nV 

u 4 N(3 3 

is implied by (3.26), (3.5) and (3.13). Finally, using (3.26), (3.7) and (3.14) we find that the 
entropy of the ideal gas of tachyons can be written as 

S t = kN ( In + In - + 4) . (3.29) 

The dependence of the entropy of a relativistic ideal gas on the temperature is shown in 
Fig. 4. We remark that, as in the nonrelativistic case, the entropies (3.27), (3.28) and (3.29) 
approach minus infinity at the absolute zero. We recall that such behavior of the entropy 
means that the classical physics fails in the limit of very low temperatures and one should 
employ quantum statistical mechanics to calculate the entropy. 

Our purpose now is to study the Helmholtz free energy. The free energy F is related to 
the partition function Z by 

F = ~hiZ. (3.30) 

Using (3.30), (3.3), (3.5), (3.7) and the Stirling formula we obtain the following formulas on 
the free energy for bradyons, luxons and tachyons, respectively: 

p \ u N Uopm J 

As an immediate consequence of the above relations, and the formula on the chemical 
potential such that 



T,V 



\dNj, 
14 




FIG. 4: A semilogarithmic graph in which entropy per particle of the ideal gas of bradyons (dotted 
line), luxons (dash line) and tachyons (solid line) given by equations (3.27), (3.28) and (3.29), 
respectively, where no = 1, k = 1 and m = 1, is plotted against the temperature. 



we get 



1 / AirVm 3 %H \ 

* = ~p r ^at + ln ) ' (3 - 35) 

" = ~p ln 4W^ (3 - 36) 

If ArcVm 3 S , 2 (u f3m) \ 



15 



Now, making use of the formula on the pressure p such that 



/' = -[%) , (-38) 

T 



dV 



and the above formulas on the free energy we arrive at the equation of state which is the 
same for bradyons, luxons and tachyons, namely 

pV = kTN. (3.39) 



Bearing in mind the formula (3.39) one may conclude, as for example Mrowczyhski [181 ]. 
that ". . . all properties of the classical gas of tachyons and the gas of bradyons are similar 
and no new phenomena have been found in the case of tachyons". We do not share such 
opinion. First of all it seems that the better version of the equation of state than (3.39) 
which does not distinguish bradyons, luxons and tachyons is the equation of state envolving 
the density of energy U/V instead of the density of particles N/V used in (3.39). Indeed, 
using (3.39), (3.10), (3.13) and (3.14) we get 

p= * / r v ( 3 - 4 °) 

where p = U /V is the density of energy and the functions /&, fi and ft corresponding to the 
case of bradyons, luxons and tachyons, respectively, are given by 

f b {u f3m) = 3 + ——— — , 3.41 

K 2 {u /3m) 

MuoPm) = 3, (3.42) 

ft{u (3m) = 3 . 3.43 

S ,2{u f3m) 

The plot of the above functions is given in Fig. 5. We point out that in the high-temperature 
limit we have (see (3.12)) 

lim f b>t {u Pm) = fi(u /3m) = 3, (3.44) 

(uo{3m)— >0 

and in the limit T — > 

fb(uo(3m) = uoj3m, uo(3m ^> 1, (3.45) 

and 

ft(u f3m) = 2, u (3m > 1, (3.46) 

following directly from (3.11) and (3.16), respectively. Furthermore, in opposition to 

16 




FIG. 5: The plot of the function fj, (dotted line), fi (dash line) and ft (solid line) given by 
(3.37), (3.38) and (3.39), related to the equation of state (3.36) for bradyons, luxons and tachyons, 
respectively, where uq = 1, and m = 1, in (3.37) and (3.39), versus kT. 

Mrowczyriski, we treat the occurrence of the maximum of specific heat for tachyons as 
their important property distinguishing them from both bradyons and luxons. Even more 
evident difference of behavior of tachyons and bradyons is the dependence of velocity on 
the temperature. Indeed, it is clear in view of the fact that upon losing energy tachyon 
accelerate, that the velocity of a tachyon should be decreasing function of temperature. Of 
course, it is not the case for bradyons. To be more specific, consider the average squared 
velocity (v 2 ). For the sake of simplicity we restrict to the preferred frame, i.e. we set uq = 1 
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FIG. 6: A semilogarithmic graph in which squared velocity of particles of the ideal gas of bradyons 
given by (3.47), where m = 1, is plotted against kT on a logarithmic scale. 

and u = 0. The average squared velocity for bradyons is given by 

V f d 3 pv 2 exp (-fiJp 2 + m 2 ^ V J d 3 p exp (-fiJp 2 + m 2 j 



2\PF 



7PF 



yPF 



mP f™x-\x 2 - l)l e - m P x dx 



, (3.47) 



where Zf b F = =1 (see (3.3)). The plot of (v 2 )^ F versus f3 is shown in Fig. 6. The average 
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FIG. 7: Average squared velocity of particles of the ideal gas of tachyons as a function of kT (see 
3.48) on a logarithmic scale. The parameter m is the same as in Fig. 6. 

squared velocity for tachyons in the preferred frame can be written as (see Fig. 7) 
V J d 3 pv 2 exp (-fiJp 2 - ™ 2 ) V Jd 



& P wfcr exp ( -f3, 1 p 2 - rri- 



mfix 



dx 



So,2(m(3) 



-, (3.48) 



where Z\f = Z ltluQ=1 (see (3.7)). 
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IV. QUANTUM RELATIVISTIC IDEAL GAS 



Our purpose now is to study the quantum ideal gas in the absolute synchronization. The 
classical formula (3.1) and the nonrelativistic relations in the case of the canonical ensamble 
[20I indicate the following form of the partition function for the relativistic gas of N particles 



z N = J2 e M-u 2 M), (4.1) 



where p° n is the energy of the system in the state labelled by n. We recall that the grand 
canonical partition function is given by 

00 

2 = Y, (4-2) 

N=l 

where s = is the fugacity and \i is the chemical potential. We point out that the product 
j3fi, and thus the fugacity, is Lorentz invariant (see (3.35), (3.36) and (3.37)). So the grand 
canonical partition function S is Lorentz invariant as well. As in the case of the classical 
ideal gas the average energy is related to the partition function by 

tf=-£l»S. (4.3) 

Furthermore, we have also 

§ = lnE. (4.4) 

From (4.4) one can obtain the equation of the state by eliminating the parameter s with the 
help of the relation 

d 

N = s— InS, (4.5) 

OS 

where N is the average number of particles of the gas. Now taking into account the statistics 
of the particles and the form of (4.2) we obtain the following form of the grand canonical 
partition function: 

3 = H(l =F sexp[-u 2 Pp°(pW\ (4.6) 

E 

where upper and lower sign refer to Bose-Einstein and Fermi-Dirac gases, respectively, and 
p°(p) is given by (2.13), (2.14) and (2.15). Using (4.6) we arrive at the following form of 
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the relations (4.3), (4.4) and (4.5): 



g = T ^ln(l T se W [-ulPp°(p)}, (4.8) 
p 

The above entities are well-defined for s > in the case of fermions and < s < 1 for 
bosons. In the limit V — > oo the sums (4.7), (4.8) and (4.9) change to integrals via the 
replacement [17J 

where we set h = 1. Hence, using (2.15), (2.16) and (2.17) we get 

U\ m 4 f°° cosh 2 tsinh 2 t ^ (4 11a) 
VJ b 2tt 2 Jq s^ 1 exp(w /3m cosht) =)= 1 

1 — [3K 2 (nu Pm) +nu pmK 1 (nu pm)]s n , (4.11b) 



2tt 2 u 2 (3 2 ^ n 

u n=l 



_ tt /" ^ xy / u 2 x 2 _|_ j / ^i„i (412a) 

V^A 47r 2 7_ 1 7 '-' a -i exp(« /W^ 2 + %|) =F 1 

- 3 e^T^ " 



tt 2 m 4 /3 4 ^ n 

u n=l 



[A m 4 f 00 cosh 2 tsinh 2 t , 

— = — / — — — dt (4.13a) 

V J t 27T s i exp(M pmsinht) =f 1 

m 2 (±l) n+1 

9 2 2R2 1J 2 — [3S 0y2 {nu f3m) - nwo/3mS , _i i i(nw /?m)]s n . (4.13b) 

u n=l 



We recall that the indices b, I, and t refer to bradyons, luxons and tachyons, respectively. In 
the case of the Bose- Einstein gases we set s < 1. The power series expansions (4.11b) and 
(4.13b) were obtained with the help of the basic properties of the Bessel functions and the 
Lommel functions presented in the Appendix. The expansion (4.12b) was derived with the 



use of the identity 21] 



o 



e bx _ c c ]ja n 
n=l 
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where T(x) is the gamma function. Furthermore, applying the same technique as with (4.11), 
(4.12) and (4.13) we find 



(P_\ 



m 



, , Ttt^ — / dtsinh t cosh t ln[l =F s exp(—u (3m cosh £)] (4.15a) 
\kT J b 2lT Uq 



m 2 r±r 



2n 2 u 2 f3 ^ 
0M n=i 



n+l 



n- 



-K 2 (nu (3m)s n , 



(4.15b) 



— ^ 

kTJi 



uq(3 
12vr 2 



dxy/ W.r 



2™2 



-i 



d|p| 



\P\ 



s^ 1 exp(u p^/u 2 x 2 + 1 \p\) =F 1 



:±i 



in+l 



71=1 



(4.16a) 
(4.16b) 



f-1 



27T 2 Mo 



/•00 

/ (it cosh 2 1 sinh i ln[l =f s exp(— n^ftm sinh £)] 
Jo 



Finally, we have 



m 



2ti 2 u 2 P 



^2 — ~ — S 0j2 (nu f3m)s r 



n=l 



m 



sinh tcosht 



2tt 2 u Jo 

,2 



m 



2ix 2 u 2 f3 



oo 

E 

71=1 



s 1 exp(M /?mcosht) =)= 1 

(±l)n+l 



eft 



n 



-K2(nu f3m)s n , 



(4.17a) 
(4.17b) 



(4.18a) 
(4.18b) 



4vr 2 



i 

dx 
1 Jo 



d|p| 



s 1 exp(M /3vn 2 x 2 +T |p|) =F 1 



1 



7T 2 ^/? 3 



y(±lT s 7 

n=l 



(4.19a) 
(4.19b) 



m 



cosh 2 1 sinh 2 1 



2ti 2 Uq Jq s 1 exp(M /5msinht) =f 1 



dt 



m 2 (±r 



E 



'1+1 



2tt 2 u 2 /3 ^-f n 

u n=l 



-S 0j2 (nu /3m)s n . 



(4.20a) 
(4.20b) 



In the preferred frame, when uo = 1, the power series expansions for bradyons (4.11b), 
(4.15b) and (4.18b) reduce to the relations originally obtained by Glaser 22|. The expres- 
sions for tachyons (4.13b), (4.17b) and (4.20b) in the particular case of the preferred frame 
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coincide up to a multiplicative normalization constant with the formulas originally derived 



by Mrowczyriski 



231 ] . We point out that an immediate consequence of (4.12b) and (4.16b) 



is the following relation 

P=\p, (4-21) 
where p = U/V. Therefore the equation of state for the quantum ideal gas of massless 
particles has the same form as in the classical case described by (3.40) and (3.42). 



A. Classical limit 

We now discuss the classical limit when N/V — > and/or T — > oo that is (3 — > 0. In view 
of (4.18a), (4.19a) and (4.20a) it is clear that in this limit the fugacity s is small as well. 
Therefore, we can approximate the series in the formulae (4.11)-(4.20) by the first terms 
linear in s. We find 

U\ m 2 



Vj b 2^uiP 
/ p \ m? 

N\ m 2 



[3K 2 (u (3m) +u (3mK 1 (u (3m)}s, s<l (4.22) 



K 2 {u f3m)s, s < 1 (4.23) 

V Jh 2^1^^ S<<1 - (424) 
From (4.23) and (4.24) we obtain the classical equation of state (see (3.39)) 

pV = kTN, s < 1. (4.25) 



We also have the relation following directly from (4.22) and (4.24) 



o , a Ki(u /3m) 
3 + Uopm- 



s < 1, (4.26) 



K 2 (u (3m)\ ' 

which is equivalent to the classical formula (3.10). Analogously, we have for luxons 

7 A 3 

s<l (4.27) 

s < 1 (4.28) 

s < 1. (4.29) 



= 1 ■ 

N\ 1 



VJ l n 2 ul(3*' 

As with bradyons we get from (4.28) and (4.29) the classical equation of state (4.25). On 
the other hand, (4.27) and (4.29) imply 

U\ 3 



nj, p< s<<1 - (4 - 30) 
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We have thus obtained the classical expression for the energy of the luxon ideal gas (3.13). 
Finally, consider the case of tachyons. The corresponding approximations can be written as 



U 
V 
(P_\ 



m 



2n 2 u 2 (3 2 



m 



[3So :2 (u f3m) - u f3mS-i : i(u f3m)]s, 



s < 1 



\kTJ t 2^^^ 



Try 



s < 1 

s < 1. 



(4.31) 
(4.32) 
(4.33) 



As in the case of bradyons and luxons we obtain from (4.32) and (4.33) the classical equation 
of state (4.25). Furthermore, using (4.31) and (4.33) we get 



1 



3 — u /3m 



S- lt i(u f3m) 
S 0i2 (u Pm) 



s < 1, 



(4.34) 



coinciding with the classical expression (3.14). 



B. The limit T -> 



Bearing in mind the technical complexity of calculations in the case of the degenerate 
Fermi-Dirac gas as well as the fact that only bosonic tachyons are admitted as a candidate 
for a dark matter we restrict to the Bose gas. Furthermore, in view of observations of [jj] 
the limit T — > in the case of Bose gas of bradyons is in fact the non- relativist ic limit and 
will not be discussed herein. Consider the ideal gas of luxons. As in the special case of the 
photon gas which is well known to be degenerate for all temperatures, we set s — 1. From 
(4.12), (4.16) and (4.19) we get 



— - rr 4 

-CT\ 
3 



P 
N 



aT 3 V, 



(4.35) 

(4.36) 
(4.37) 



where C = 3((4)k 4 /7i 2 u 4 = n 2 k 4 /30u 4 , and a = ((3)k 3 /n 2 u 4 ; ((x) = J2n=i & is the 
Riemann zeta function. In the case of the photon gas and the preferred frame (u = 1), the 
constant C and the constant a, taking into account the two polarization states of a photon, 
should be multiplied by 2, i.e. C is then the Stefan-Boltzmann constant and a = a p hoton, 
respectively. 
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Now the entropy can be written as (see (3.24)) 



S = k\n~ + ^. (4.38) 

The equations (4.38), (4.4), (4.35) and (4.36) taken together yield 

S = ^VCT 3 . (4.39) 

As expected, the entropy vanishes at zero temperature. The formula (4.39) can be also 
obtained from the well-known relation 

S= [ T ^dT, (4.40) 



o 



where 



following directly from (4.35). 

We now study the Bose gas of tachyons in the low-temperature limit. In this limit one 
cannot simply replace sums (4.7), (4.8) and (4.9) by integrals because when s ~ 1 the terms 
referring to the ground state can give finite contributions to series. Therefore, we separate 
the sum (4.9) into two contributions — the number of particles in the ground state Nq and 
the number of particles in the excited states 20]. Using (4.20b) we find 

n 



2 °° -i 

N = N ° + V lTTT-a E So,2(nuoPm)s n , (4.42) 

u n=l 



where 

N = (4.43) 
s 1 — 1 

We recall that the accumulation of bosons in the ground states is known as Bose-Einstein 
condensation. Taking the limit of the strong degeneration s = 1 and utilizing for T — > 0, 
that is /3 — > oo, the asymptotic formula (A. 10) we obtain from (4.42) the relation 

From (4.44) it follows that the Bose-Einstein condensation occurs at temperatures lower 
than the critical one given by 

Wt^' (4 - 45) 
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and densities higher than the critical density such that 

mk 2 



d. 



12m 3 



T 



(4.46) 



Furthermore, using (4.44) and (4.45) we get 



No 
N 



T < T r . 



(4.47) 



For temperatures greater than T c but close enough to T c to enable setting s = 1, we have 

U 



m 



V " " 7r 2 M ^ 3 " C(3) ' 
p m 



kT 2n 2 u\^ 2 
N m 



C(3), 
C(2). 



(4.48) 
(4.49) 

(4.50) 



V ' 2ir 2 u 3 (3 2 

Hence, taking into account (4.49) and (4.50) we arrive at the equation of state of the form 



(4.51) 



Using (4.48) and (4.50) we also obtain 



Therefore the specific heat is 



U = Z^NkT, T > T c . 



C v = 2^|| Nk, T > T c . 



(4.52) 



(4.53) 



For temperatures lower than T c one should use N — N instead of N. Therefore the equation 
of state takes the form 



and (4.52) is replaced by 



P 



U 



m 



C(3)fc 3 T 3 , T < T c , 



2tt 2 



ii'r 



(4.54) 



m 



C(3)rT 3 \/, T < T c . 



(4.55) 



o 



Using (4.38), (4.4), (4.54) and (4.55) we obtain the following formula on the entropy: 

S = -^UC(3)A; 3 T 2 , T < T c , 



(4.56) 
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implying vanishing of the entropy at zero temperature. The relation (4.56) is also a conse- 
quence of (4.40) and the expression on the specific heat such that 

tTT) 

C v = V((3)k 3 T 2 , T < T e , (4.57) 

following immediately from (4.55). The discussion of the relations satisfied by the thermo- 
dynamic functions for the non-relativistic Bose gas in the case of T > T c , and T < T c can 
be found in the book 24|. 



V. CONCLUSIONS 

In this work we have derived the Lorentz covariant form of thermodynamic functions for 
the relativistic ideal gas in both classical and quantum cases. We stress that the applied 
approach based on the concept of the preferred frame is the only one which enables for- 
mulation of covariant statistical mechanics of tachyons. On the other hand, an advantage 
of the formalism introduced in this paper is that it allows to study from the unique point 
of view all kinds of relativistic gases including bradyon, luxon, and tachyon ones. Bearing 
in mind the possible applications of the observations of this work we point out the specific 
heat of the tachyon gas showing behavior analogous to Schottky anomaly (see Fig. 3), and 
the formula (4.43) on the critical temperature for the Bose condensation which seem to be 
of importance for the determination of the mass of a tachyon. The results obtained in this 
paper can be also helpful in discussion of the dark matter and dark energy in astrophysics 
and cosmology. Indeed, the observations of stars motion in galactics, galactic clusters, clus- 
ter masses as inferred from gravitational lensing strongly suggest existence of an exotic dark 
matter component of the universe. Moreover, from observations of supernova IA popula- 
tions we know that the universe accelerates which demands existence of exotic dark energy. 
A possible attempts to explanation of the exotic content of the universe need some radical 



extension of standard physics. One of candidates is tachyonic fluid (see for example 26]) 
or the rolling tachyon field 0, \^\ . 
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APPENDIX 

We first recall some properties of the modified Bessel functions (Macdonald functions) 
K v (x). These functions have the integral representation of the form 

K v (x) = / e~ xcosht cosh vtdt, x>0. (A.l) 
Jo 

They satisfy the following recurrence relations: 

K'„{x) = V -K v (x) - K v+X {x) = ~K v {x) - K v _i(x) } (A.2) 
where prime designates differentiation. We have the asymptotic formulas 

K n {x) = i( n -l)!0) , (A.3) 

K v {x) = J^e~ x , x » 1. (A.4) 



We now briefly sketch the basic properties of the Lommel functions S tl}U (x) 2l|, |25j. The 
Lommel functions Sq )V (x) are given by 



v x roc 



S 0tU (x)= e~ xsmht cosh vtdt = - / e~ xsmht coshtsmh vtdt, x > 0. (A.5) 
Jo v Jo 

The Lommel functions S UyU (x) can be expressed in terms of the Struve functions H u (x) and 

the Bessel functions Y u (x) (Neumann functions) also designated by N u (x). We have 

S VtV (x) = 2^ 1 v ^i> + §)[H„(x) - Y v (x)\. (A.6) 

The recurrence relations for S ll ^ v {x) are of the form 

S^vOc) = - [fa + l) 2 - u 2 ]S^{x), (A.7) 
S'w( x ) = ~ S ^A X ) + (L L -v- 1)5^-1,^+1 (x) = ~S^ u (x) + (p + v- l)S,i-i, v -i(x). 

(A.8) 

For small and large argument So^{x) can be approximated as 

SoM = ^ x < 1, (A.9) 
S , 2 (x) = -, x > 1. (A.10) 

3j 
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We remark that So,2(x) and K 2 (x) are approximated by the same function in the limit 
x < 1. 
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